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Abstract 

We study semiclassical approximations to the time evolution of coherent states 
for general spin-orbit coupling problems in two different semiclassical scenarios: The 
limit h ^ is first taken with fixed spin quantum number s and then with hs held 
constant. In these two cases different classical spin-orbit dynamics emerge. We prove 
that a coherent state propagated with a suitable classical dynamics approximates the 
quantum time evolution up to an error of size \^ and identify an Ehrenfest time scale. 
Subsequently an improvement of the semiclassical error to an arbitray order h^^'^ is 
achieved by a suitable deformation of the state that is propagated classically. 
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1 Introduction 



Ever since their introduction by Schrodinger as early as 1926 |Sch26j . coherent states have 
found an increasing range of apphcations in quantum mechanics, see e.g. jKS85[lPer86j . In 
a semiclassical context their virtues become particularly transparent in attempts to relate 
the quantum time evolution of a system to its classical trajectories. Coherent states can, 
e.g., even be used to identify the limiting classical dynamics of a given quantum system. 

However, apart from the exceptional case of the harmonic oscillator that Schrodinger 
chose for his construction, every quantum wave packet necessarily disperses. Schrodinger's 
original intention to mimic classical trajectories in quantum mechanics can therefore only 
be put into practice up to the time scale on which wave packets begin to delocalise. Beyond 
that the quantum time evolution looses its tight relation to classical trajectories, although 
coarser classical structures possibly remain to be of influence [SBO?, SchMj. 

More recently the notion of an Ehrenfest time was introduced ^Chi79|IZas81j . intended 
to indicate that the Ehrenfest relations can only connect quantum dynamics and classical 
trajectories on limited time scales. For classical dynamics with positive Lyapunov expo- 
nents it is argued that the Ehrenfest time is logarithmic in h. This conclusion can be drawn 
from the observation that coherent states are localised in phase space on a scale of \/h, and 
that an unstable classical dynamics expands domains in phase space with exponential rates 
in the unstable directions. Thus for times beyond 2^ | log h\ a coherent state is no longer 
localised in directions that are expanded with an exponent A. A finer analysis reveals that 
the precise value of the Ehrenfest time depends on the problem that is studied; e.g., using 
L^-norms to measure the difference between the quantum time evolution of a coherent state 
and a coherent state that is propagated with the classical dynamics, a critical time scale 
of g^l log^l was proven to hold |.CR97j . On the other hand, the same difference measured 
in terms of expectation values of observables can be controlled up to times of the order of 
^1 log^|. For details see |CR97| IBBOOt IBR02j . It can moreover be shown that on finite 



time intervals a coherent state is exponentially localised around the corresponding classical 
trajectory |H.inn]. 

Except for heat kernel asymptotics in the case of particles in non-abelian gauge fields 
|HPS 83j most of the previous work on a semiclassical control of the propagation of coherent 
states is concerned with systems that possess only translational degrees of freedom. In 
this article it is our aim to extend these investigations to systems with non-relativistic 
spin-orbit interactions. After having identified appropriate coherent states, we intend to 
compare solutions of the Schrodinger equation 

dip 

= Hij{t,x) , 

where the initial wave function ?/'(0) is a coherent state, with a coherent state that is 
evolved along suitable classical trajectories. The quantum Hamiltonians that we wish to 
allow are of a general spin-orbit coupling type, 

H = Ho{Q,P) + CiQ,P)-S , (1.1) 
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with Q, P, and S denoting the standard position, momentum, and spin operators, re- 
spectively. Examples of such Hamiltonians arise when the spin is coupled to an external 
magnetic field, such that G = ^-B, or in the context of atomic spin-orbit coupling with 
C being proportional to orbital angular momentum. 

Apart from atomic and molecular physics spin-orbit coupling also plays an important 
role in nuclei, where it essentially determines their shell structure (BMGOj, as well as in 
solid state physics. In the latter case recent experimental progress towards controlling the 
spin dynamics of electrons in semiconductors SFHZOl calls for a theoretical description of 
such set-ups. As opposed to some pure quantum calculations semiclassical considerations 
are often particularly transparent and provide a clear physical picture. With our work 
we therefore intend to improve the understanding of spin-orbit coupling by establishing 
mathematically rigorous statements about the quantum dynamics of localised particles 
with spin and their relation to appropriate classical trajectories. 

One issue to be settled is how the semiclassical limit should be performed in the pres- 
ence of spin-orbit interactions. In principle two parameters controlling the passage to a 
classical description are available, which are associated with the two types of degrees of 
freedom: translational and spin. On the one hand, with (an effective) h approaching zero 
the semiclassical limit is achieved in a standard way for the translational degrees of free- 
dom. On the other hand, for an isolated spin h can be eliminated from both kinematics 
and dynamics. The role of a semiclassical parameter is then taken over by 1/s, where 
s = 1/2,1,3/2,... denotes the spin quantum number. When both types of degrees of 
freedom interact through a spin-orbit coupling one can therefore pass to the semiclassical 
limit in various ways. In the absence of a theory that is uniform in both h and 1/s we 
subsequently focus on two important scenarios: 

The most straight forward approach is to view, say, an electron as a particle with fixed 
spin 1/2 and to employ h as the only semiclassical parameter. In the limit ^ — >■ the energy 
scale of the translational part Hq in (jl.ip then dominates that of the spin-orbit coupling 
term, since in the latter the spin operator S is proportional to h. Although it might appear 
that thus the spin has evaded the leading order semiclassical description, it does in fact 
contribute in an essential way through a classical spin precession driven by the orbital 
motion, see |BK99a[ IBC^OO^ IBC^KOH IB(;n4j . E.g., in classically chaotic systems this type 
of spin motion is responsible for the quantum eigenvalue spectrum to possess correlations 
of the Gaussian symplectic ensemble of random matrix theory |BK99bj . Moreover, in this 
semiclassical framework the exact spectrum of the relativistic hydrogen atom is recovered 
|Kep03| , and anomalous magneto-oscillations in semiconductor devices can be described 
to a good approximation |KWn2j . 

A second option for the semiclassical limit is to keep the "classical spin" hs at a fixed 
value S, thus performing h and s ^ oo simultaneously. In this scenario the energy 
scale of spin-orbit interactions remains comparable to that of the purely translational part, 
leading to a classical spin-orbit Hamiltonian. Therefore, coupled Hamiltonian dynamics 
emerge with classical particle trajectories infiuenced by the spin. This scenario enables an 
immediate classical description of the Stern- Gerlach experiment, and generally corresponds 
to a "strong" spin-orbit coupling. 
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In this paper we examine the propagation of coherent states under the influence of 
spin-orbit interactions in both of the above mentioned semiclassical scenarios. In section |21 
we first provide a precise characterisation of the quantum Hamiltonians under investigation 
and then describe the classical dynamics that will result in due course. Section IHl is devoted 
to outlining the construction of coherent states for both translational and spin degrees of 
freedom, along with their basic properties. Our principal results are developed in section |3J 
For both semiclassical scenarios separately we extend the approach devised previously 
[HclTSl iLitSGj, ICR97j in systems without spin in that we first construct suitable approximate 
Hamiltonians that propagate coherent states exactly along classical trajectories. We then 
prove that, measured in Hilbert space norm, the full quantum dynamics differs from a 
classically propagated coherent states by an error of size ^Jfi as long as finite times are 
taken into account. The vanishing of this difference up to some, semiclassically infinite, 
Ehrenfest time is also established. Subsequently we improve the semiclassical error to 
0{h^^'^) for arbitrary G IM by replacing the classically propagated coherent states with 
a suitable sum of squeezed states. Again such a procedure is possible up to the Ehrenfest 
time. We conclude in section El with discussing some implications of our main results. 

2 Background 

It is our aim to investigate the time evolution of an initial coherent state in both transla- 
tional and spin degrees of freedom generated by a general spin-orbit quantum Hamiltonian 
with an emphasis on a semiclassical description. This is the reason why we represent 
quantum observables as matrix valued semiclassical pseudodifferential operators within 
the framework of Weyl calculus, see |R,ob87[ IDS99j for details. The quantum Hamiltonians 
H under consideration are defined on the domain ^(R"') ® C^'^^^ in the Hilbert space 
L'^iW^) ® C^''^^ and are of the form 

{Hi,){x) = (op^[i7]^)(x) = j j^^^^^m^-y) H[^,i)^{y) dyd^ . (2.1) 

Here T*E'^ = E*^ x E'^ denotes the cotangent bundle over the euclidean configuration space 
E'^, i.e. the phase space of the translational degrees of freedom. The spin s = 1/2, 1, 3/2, . . . 
is described by the matrix degrees of freedom of the Weyl symbol H and will later be 
represented on its phase space S^. Spin-orbit Hamiltonians are characterised by symbols 
of the form 

H{x, = H^{x, + h Cix, ■ dnsicT/2) , (2.2) 

where Hq and the components Ck, k = 1,2, 3, of C are real valued and smooth functions 
on T*R'^ which for all multi- indices a and f3 satisfy the growth estimate 

\d:d',F{x,0\ < K^, (1 + \x\'Y"^'\l + (2.3) 
with suitable constants > and M = (M^, M^) G R^. 
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The spin-orbit couphng term in ()2.2|) contains the spin operators 

Sk ■■= hdTTs{crk/2) , k = 1,2,3 , 
obeying the well known commutation relations [Sk, Si] = ihetim Sm- Here 

are the Pauli matrices, considered as elements of the real Lie algebra su(2), and dTr^ denotes 
the (2s+l)-dimensional representation of su(2) derived from the corresponding unitary irre- 
ducible representation tt^ of the Lie group SU(2) according to A'Ks{X) = i^ 7rs(e"'^^) 

The time evolution U (t) = e~'^^* generated by the quantum Hamiltonian will be unitary 
provided that H itself is essentially self-adjoint on the domain C^(M!^) ® C^*"*"^. In the 
present framework this is guaranteed, for sufficiently small h, once the symbol H is such 
that if + i is elliptic, i.e. if 

\\{H{x,0 + ir'\\<c{l + \x\y^^\l + \^\y'^' (2.5) 



holds for all (x,^) G T*R'^ with some constant c > and M as in ()2.3jl : here || ■ || is an 
arbitrary matrix norm. Details can be found in |Rob87irn^ . In the following we assume 
this condition to hold and do not notationally distinguish between H and its self-adjoint 
extension. 

In the semiclassical limit we will have to deal with two types of classical spin-orbit 
dynamics. In the first case only the translational degrees of freedom evolve under a Hamil- 
tonian fiow. This is defined on the phase space T*R'^ and is generated by the classical 
Hamiltonian Hq. Thus ^Q{q,p) = [q(t),p(t)) satisfies Hamilton's equations of motion, 

q{t) = d^Ho{q{t),p{t)) and pit) = -d,Ho{qit),pit)) , 

with initial conditions (g(0),p(0)) = {q,p)- This fiow then drives a classical spin through 
the equations of motion 

h{t) = C{qit),p{t)) xn{t) 

on the sphere with initial condition n(0) = n. Here n G R'^ with |n| = 1 is considered 
as a point on S^. The curve n{t) therefore describes the Thomas precession of a normalised 
classical spin vector on along the trajectory in T*E''. The combined dynamics 

{q,p,n) ^ [¥Q{q,p),n{t-q,p,n)) (2.6) 

yield a fiow on the product phase space T*R'^ x S^, which is a symplectic manifold whose 
symplectic form is composed of the natural symplectic forms of its factors. This fiow has 
the form of a skew product, see |(yFS82j for details, and thus is not Hamiltonian; however, 
it leaves the natural volume measure derived from the symplectic form invariant. 

The second fiow relevant for our subsequent discussion includes a classical spin dynamics 
coupled to the motion of the translational part in a Hamiltonian manner and is also defined 
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on the product phase space T*E'^ x S^. These dynamics are generated by the classical spin- 
orbit Hamiltonian 

HUx, e, n) := Ho{x, + Sn ■ C{x, , (2-7) 

where the constant S > measures the length of the classical spin vector s := Sn. The 
Hamiltonian flow ^l^(q,p,n) = n(t)) is therefore determined by the equations 

of motion 

q{t)=diiH,o{q{t),p{t),n{t)) , 

p{t) = ~d,H,o{q{t),p{t),n{t)) , (2.8) 
h{t) = C{q{t),p{t)) X n{t) . 

The Hamiltonian coupling of the degrees of freedom prescribed by these equations imply 
that in contrast to the previous case the translational dynamics are affected by the spin. 

Apart from the associated classical flow in semiclassical approximations of quantum 
dynamics also the linear stability of the flow plays a role. Quantitatively this can be 
measured in terms of the Lyapunov exponents, see the appendix for a discussion. They 
express the rate of phase space expansion or contraction, respectively, induced by the flow 
in different tangent directions. Moreover, the differential of the flow is a symplectic map on 
the tangent bundle of phase space. Its metaplectic representation is an essential ingredient 
in the semiclassical propagation of coherent states. 



3 Coherent states 

Within the setting outlined in the preceding section we wish to describe the time evolution 
of an initial coherent state semiclassically. The starting point therefore is the Schrodinger 
equation 

dip ^ 

i^-^(i,a;) = Hi){t,x) with ^/'(0,x) = (v?J^p) 'S)(pn){x) , 

whose initial condition is the product of a coherent state V^^p) of the translational degrees of 
freedom and a spin-coherent state 0„. The principal question we then address is to what 
extent the quantum mechanical time evolution can be approximated by some classical 
dynamics, i.e. we want to estimate the difference 

in terms of h, where (g(t),p(t), n(t)) is an appropriate classical trajectory and e'"*-*^ is a 
suitable phase factor. 

For both types of coherent states, ffgp-j G L'^(R'^) and 0„ G C^'^"*""'^, we use Perelomov's 
construction |Per86j that applies to a general Lie group G with unitary irreducible represen- 
tation TT on a Hilbert space 7i: Fix a non-zero vector ^ ^ consider := 7i{g)'$Q for 
every g E G. Hence the vectors \l'p and deflne the same quantum state, i.e. = e^°"$g, 
if and only if g~^h lies in the stability subgroup H G G of the vector \E'o) 

H:={geG; n{g)^o = e^'^'^^o} • (3.2) 
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The quantum states generated by the vectors '^g, g E G, can thus be labeled by the points 
T] of the coset space G/H. A section g{ri) in the bundle G ^ G / H then determines a 
choice of vectors 

% := = 7r(r/(r/))^o , for r/ G G/H , (3.3) 

representing these states. The vectors $^ are called coherent state vectors for {G,tt,H). 

The two types of coherent states that play a role in the present setting can be con- 
structed according to this general scheme by choosing the Heisenberg group G = if(R'^) 
for the translational part and the group G = SU(2) for the spin part. We now describe 
the two situations that emerge from this procedure separately. 



3.1 Coherent states for the Heisenberg group 

The Heisenberg group if(R'^) is a non-compact (2(i+l)-dimensional Lie group that consists 
of the elements {q,p, A) with {q,p) G T*E'^ and A G IR. The group multiplication is given 
by 

{q,p,X) {q',p',X') = {q + q',p + p',X + X' + l{pq' -qp')) . 

According to the Stone-von Neumann Theorem any unitary irreducible representation vr of 
HiW^) that fulfills 7r(0, 0, A) = en^ is unitarily equivalent to the Schrodinger representation 
pn on L2(E^), 

(pft(g, p, X)4j) (x) = e^^ (e^(P<3-gP)^) (3.) = gi(A+p(x-|q)) _ _ 

Here Qk and Pk, k = 1, . . . ,d, are the standard self-adjoint position and momentum oper- 
ators defined on suitable domains in L^(E'^). 

In order to construct coherent states for the Heisenberg group we therefore consider the 
Schrodinger representation pfi on L^(R°'). One immediately sees that given any non-zero 
vector \E^o ^ its stability subgroup is if = {(0, 0, A); AG E}. Thus coherent states 

can be labeled by points {q,p) G G/H, i.e. by points in the phase space T*E'^. This labeling 
can be achieved in terms of the section g{q,p) := {q,p, —^qp) in G ^ G/H. One usually 
prefers a reference vector \E'o ^ ^^^(IR'^) that is normalised, rapidly decreasing, and satisfies 

(*o,Q^o) = and (^o,P^o)=0, 

so that any reasonable lift of this vector to the phase space T*R'^ is concentrated at 
(0, 0) G T*R'^. A convenient choice with these properties is 

^o^(^) -(^(detImi?)V^e-^^% 

where B is some complex symmetric d x d matrix with positive-definite imaginary part. 
The coherent states ()3.3|1 that follow from the above definitions now read 

vfq,p)i^) = {Ph{q,P,-lqp)^Q){x) 

= (detIm5)i/^eK^'(-'?)+^(-^)^(--'^)) . ^^"^^ 
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Note that these coherent states differ slightly from more conventional choices for which 
B = ilrf and the section g{q,p) = {q,p, 0) are used, leading to a different phase convention. 
Despite the fact that after allowing for more general matrices B the coherent states loose 
the minimum uncertainty property, this generalisation will prove useful since the action 
of the metaplectic representation on them can be conveniently expressed in terms of B, 
see also |Sch01j . The alternative phase convention is of less consequence but simplifies the 
notation. 

Although from the above construction it is obvious that a coherent state V'^p) is con- 
centrated in some neighbourhood of the point {q,p) in phase space it is instructive to 
calculate explicit phase-space lifts. E.g., its Wigner transform is given by 



Jh'^ [o.o) 

^ 2d^-^{{x,0-{q,p))GB{{x,0-{q,p)) ^ 

where is the positive-definite symmetric 2d x 2d matrix 



Gn :-- 



ImB + ReB{lmB)-^ReB -ReB{lmB)-^ 



^' » -{ImBy^ReB (ImB)-^ 

This representation reveals a concentration of the coherent state in the vicinity of the 
phase-space point {q,p). Moreover, since the sum of position and momentum uncertainties 
reads 



the spreading of the coherent state in phase space can be measured in terms of Gb- 



(3.6) 



3.2 Spin-coherent states 

In quantum mechanics the spin of a particle is implemented through the (2s+l)-dimensional 
irreducible representation tt^ of the compact Lie group SU(2), where s = 1/2, 1,3/2, . . . 
denotes the spin quantum number. Within Perelomov's framework spin-coherent states 
are hence constructed from (SU(2), vr^, C^*"*"^). The reference vector \I'o G C^''"'"^ can be 
chosen such that the coherent states possess the minimum uncertainty property; this is 
achieved with \E'o being a maximal weight vector for the irreducible representation dvr^ of 
the Lie algebra su(2). 

The real Lie algebra su(2) consists of the hermitian and traceless 2x2 matrices X, 
such that e~^^ G SU(2). A convenient basis of su(2) is formed by the Pauli matrices ()2.4j) . 
We also consider the complexified Lie algebra su(2)c := su(2) (g) C with basis given by 

X± := ^{(Ji ± i(T2) , X3 := ^0-3 , 
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and commutation relations 

[X3,X±] = ±X± , [X+,X_] =2X3 . 

The vector X3 spans a Cartan subalgebra, which exponentiates to a maximal torus T ~ 
U(l) in SU(2), and X± e su(2)c span the root spaces q± C su(2)c. Their representations 
d7rs(X-|-) are raising and lowering operators, respectively. More precisely, the representation 
space C^^"*"^ decomposes into a direct sum of the one dimensional eigenspaces of dnslX^) 
(weight spaces) = {0 G C^*^"'^; d7rs(X3)0 = m0}, where m = —s,—s + l,...,s. 
The raising and lowering operators d7rs(X±) map the weight spaces into one another, 
d'n's{X±)Vm = Vm±i for m 7^ ±s. The weights m = ±s are called maximal and minimal 
weights, respectively. The corresponding weight vectors are annihilated by the raising or 
lowering operator. In the usual angular momentum notation a normalised weight vector is 
denoted as |s, m). 

For a given representation vr^ of SU(2) we choose a maximal weight vector |s, s) as the 
reference vector \E'o- According to ()H.2|1 the stability group of this vector is 

H = {g = e-'^^^- AG [0,27r)} = U(l) , 

which can be identified with a maximal torus T. Thus coherent states are labeled by points 
in the coset space 

G/H = SU(2)/U(1) ^ . 

As in the case of the Heisenberg group this manifold is naturally symplectic and can be 
viewed as the corresponding classical phase space. 

The definition of coherent states finally requires a section in G G/H, i.e. in the Hopf 
bundle SU(2) — > S^. This principal U(l)-bundle, however, is non-trivial so that no smooth 
global section exists. We therefore here give local constructions that, nevertheless, allow 
for suitable interpretations in terms of global objects. We parameterise points on by 
n G with \n\ = 1 and use spherical coordinates, n{9,ip) = (sin ^ cos sin ^ sine/?, cos 6*) 
with 9 G [0,7r) and (p G [0,27r). Introducing := (— siny9, cosy^, 0) our choice of a local 
section reads (see also |Per86j ) 

, _^-he., <r_( cosf -sinfe-'^\ 
^"-^ - ^nfe'^ cosf ) ■ 

Under the double covering map R : SU(2) S0(3) that is defined through {R{g)x) ■ cr = 
gx ■ crg~^, the matrix gn{e,ip) corresponds to the rotation R{gn{e,Lp)) about the axis with 
angle —9, such that R{gn)e3 = n, where 63 = (0,0,1) represents the north pole on S^. 
With these choices spin-coherent states are the normalised vectors 

(Pn ■= T^s{9n)\s,s) . (3.7) 

These states are conveniently represented on the phase space through the Husimi trans- 
form (see |Per86j ). 

h[4>n]{m) := \ {4>m,(pn)\ = 



m ■ n 
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which clearly indicates a concentration, in the semi classical limit s — > oo, of 0„ at the 
point n G S^. 

Our next aim is to investigate the relation between the propagation of a spin-coherent 
state ()H.7j) generated by a (time-dependent, linear) spin-Hamiltonian 

4pin = C{t) ■ S (3.8) 

defined on C^'^"'"^, and a suitable classical time evolution n{t) on S^. Here S denotes the 
vector of spin operators Sk = hdns{ak/2). The dynamics of a coherent state (pn follows 
from the equation 

i;z^(t) = 4pin0(t) with 0(0) = 0„. (3.9) 
A solution of this problem can be related to the curve g(t), t G R, in SU(2) determined by 
g{t) + ^Cit)-(Tgit) = with ^(0)=idsu(2) (3.10) 

through 

0(t) = TTs{g{t))<f){0) = ns{g{t)gr,)\s,s) . (3.11) 

An associated classical time evolution then arises from the adjoint action of g{t) on n ■ cr G 
su(2) via n{t) ■ cr := g{t)n ■ crg{t)~^ = {R{g{t))n) ■ cr. This implies 

n{t) = C{t) X n{t) with n(0) = n . (3.12) 

The corresponding coherent state vector (f)n{t) differs from the quantum time evolution (j){t) 
of (pn only by a phase; both vectors therefore describe the same quantum state. 

Since this phase is required for later purposes, we now determine it explicitly. To this 
end we notice that n{t) can on the one hand be represented as 

n{t)=R{g{t))n = R{g{t)g^)e3 , 

and on the other hand as 

n{t) = R{gn{t))e3 ■ 

Thus, under the double covering map, g~^^^g{t)gn G SU(2) is associated with a rotation 
about 63 with some angle g{t), such that 

9;.lt)9it)gr. = e^^W-^ G T . 

From ()3.11|) it now follows that 

0W = 7r,((7(t)(7„)|s,s) = ns{gr,it))^s{e'^'^'^''')\s,s) = e^^^« , (3.13) 
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thus confirming the claimed relation between the quantum and 'classical' propagation of 
the spin-coherent state Due to the explicit dependence of the phase on s it suffices to 
calculate the angle g(t) for s = |. For this one notices that in polar coordinates 

i>{t) = ei««> ,„,„|1. i> = (^^^ffl£!l,^^ . (3.14) 

In a standard calculation (see e.g. jBK99b] ) g(t) can now be determined by using ()H.14j) in 
equation leading to 

g{t) = - ^ (^C(t') ■ n{t') + (1 - cos e{t'))if{t')'^ dt' . (3. 15) 

If one introduces a classical spin vector s := Sn, with some S > 0, one can relate the angle 
g{t) to Hamilton's principal function of the spin. The observation that 

LsMt) = -C{t) ■ s{t) - 5 (1 - cose{t))^{t) 

is the Lagrangean of the classical spin motion implies Sg{t) to be the spin-action Rspin{t). 

4 Time evolution of coherent states 

In this section we discuss the time evolution of coherent states in two different semiclassical 
limits. In the first scenario we consider h while the spin quantum number s is fixed. 
This will imply that primarily the translational degrees of freedom become semiclassical. 
The spin-orbit interaction therefore occurs on the level of the subprincipal symbol of the 
Hamiltonian ()2.2|) . enforcing the skew-product structure ()2.6p of the resulting classical 
dynamics with the translational motion driving the spin. 

In the second scenario we fix the product S := hs and hence consider the combined 
limits ^ — > and s oo. Thus both types of degrees of freedom are treated semiclassically 
on equal footing. This results in a classical spin-orbit coupling with the Hamiltonian 
dynamics ()2.8|) generated by the function ()2.7|) . 

We begin with the first scenario which is close to the time evolution of coherent states 
without spin degrees of freedom. 

4.1 Semiclassics with fixed spin 

In the present scenario h is the only semiclassical parameter so that we consider the quan- 
tum Hamiltonian ()2.1|) as a Weyl operator with matrix valued symbol ()2.2|) that has a 
scalar principal part; the subprincipal symbol then contains the spin-orbit coupling. This 
setting ensures that the propagation of coherent states is closely analogous to the case 
without spin, compare |(yR97j . 

Guided by this analogy we first construct an approximate Hamiltonian that propagates 
coherent states exactly. Regarding the translational part we exploit the fact that the time 
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evolution generated by a quadratic Hamiltonian preserves the form ffgp-j given in ()3.4|) 
of a coherent state for the Heisenberg group. The spin part of the coherent state shall 
be propagated by a Hamiltonian of the form ()3.8j) and can hence be calculated explicitly. 
Using the convenient notation := (x, ^) G T*R°', we now consider the Taylor expansion 
of the symbol ()2.2|1 about some smooth curve z{t) = in phase space. The Weyl 

quantisation of the leading terms in the Taylor expansion (of different order in the principal 
and in the subprincipal symbol), 

Hgit^w) := l^Hl;'\zit)){w - zit)y + hC{zit)).dn,icT/2) , (4.1) 

kl=o 

yields a quantum Hamiltonian Hqit) that is quadratic in Q and P and linear in S. Here 
h'q\w) stands for the derivative 5^ifo(w) of order |z/| in the 2d components of w; = (x, ^). 
The time evolution ^Jqit) G /.^(E"^) ® C^^"*"^ of a coherent state V'^p) ® <Pn generated by 
the approximate Hamiltonian, 

ih^it) = HQ{t)ijQ{t) with MO)=vL)^^r., (4.2) 
can be expressed in terms of a coherent state: 

Proposition 4.1. The solution of the quadratic Schrddinger equation \4-^ is a time- 
dependent coherent state with an additional phase, 

^Q(t) = e V ^' 2 yjj (^^^J/^^^^^^^ ® 0^,) . (4.3) 

Here (g(t),j9(t)) = ^Q{q,p) is the solution of Hamilton's equations of motion generated by 
the principal symbol Hq, 

q{t) = d^Ho{q{t),pit)) , Pit) = -d,Ho{qit),p{t)) , (4.4) 

with initial condition (g(0),p(0)) = {q,p) and principal function 

Roit) = !\p{t')q{t') - HMt').P{t'))) dt' . (4.5) 

^0 

The complex symmetric d x d matrix B(t) is given by 

B{t) = {d,q{t)B + d,q{t)) {d,p{t)B + d,p{t))-' , (4.6) 

where the derivatives are taken with respect to the initial conditions; it also gives rise to 
the Maslov phase cr{t). Moreover, n{t) is a solution of the spin precession equation i3.1^} 
in which C(t) stands for C (^q(t) , p{t)^ from \4-^}) >' Qi^) then is the associated angle 
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Proof. For the proof we adapt the method of |Sch01j to the present situation and therefore 
introduce the ansatz 



nit) 



in equation ()4.2|) . To deal with the spin contribution to the left-hand side we use the fact 
that according to ()3.9p and ()3.13p 



d 



"^^V- ^r^{t)) = hC{q{t),p{t)) ■d'Ks{(T/2)e''-- (pr^(t) , 

if and only if n{t) solves ()3.12j) . It hence remains to consider (see (ESEqI]) 



(4.7) 



dt 



Ho + H',^,{x - q{t)) + H',^^ . B{t){x - q{t)) + -{x - q{t)) ■ H'^^,,{x - q{t)) 



+ ^(x - q{t)) ■ H'l^,B{t){x - q{t)) + ^(x - q{t)) ■ B{t)Hl^,{x - q{t)) 
+\{x - qit)) . B{t)H'l^^B{t){x - q{t)) + | tr«,^, + H'l^^B{t)) 



Here the abbreviations ifg.x = dxHo and H^^^ = d^dxHo, etc. have been employed. These 
expressions are to be evaluated at z{t). Comparing coefficients of powers of h and of 
(x — q(t)) in ()4.8|) then yields the conditions 



e = qp- Ho , 
-p + Bq = H'^^^ + BH'^j^ , 

-B = Hq^^^ + H'^^^B + BH'^^^ + BH'^^^B 



(4.9) 



7 



1= t,[H'^^^^ + H'l^^B)+\sQ 



7 



With the identification Rq = 9 the first and the second equation immediately imply ()4.5|1 
and ()4.4j) . respectively. 

The other two equations involve the time evolution B{t) of the complex symmetric 
d X d matrix B with positive-definite imaginary part; they determine the action of the 
metaplectic group on the vector ffqpy At this stage we recall that the symplectic group 
Sp{d, E) acts on the Siegel upper half-space (see |Fol89j ) 

Ed := {Z G Md(C); = Z, ImZ > 0} 



via 



S[Z] = {SuZ + Si2){S2iZ + ^22)-' , where S = f ^'A G Sp(rf, E) . 

^021 022/ 
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In the present context the differential of the Hamiltonian flow $q generated by the classical 
Hamiltonian Hq is symplectic, S'o,z(t) := D^q{z) G Sp((i, E), and hence can act on the 
initial value B G S^. Indeed, 

5(t) = 5o,.(t)[5]GS, (4.10) 
yields the solution of the third equation in ()4.9|) and implies fl4.6p . The fourth equation 

— 1 /2 

requires the introduction of the Maslov multiplier m{S,Z) := (det(S'2i^ + 522)) for 
S G Sp((i, E) and Z G S^. This allows us to define the Maslov phase a(t) through 
gifo-(t) _ fn{0{t),it), where 0{t) is an orthogonal symplectic matrix that is uniquely 
associated with So,z(t). One can then show (cf. |Sch01j ) that 

7(t) = (detIm5(t))'/'e^^'^W+^'^^W. 

□ 

We remark that the state ()4.Hj) is closely analogous to the respective solution without 
spin-orbit coupling. It differs from the latter only by the factor e^'^^'^*^0„(t). This observation 
not only means that quantum mechanically the translational part and the spin part are 
not entangled, but also on the classical level the translational dynamics are independent 
of the spin precession n{t). The combination of classical translational and spin motion 
rather has the structure of a skew product ()2.(jj) . indicating that only the spin dynamics 
depends on the translational part, and not vice versa. 

Our aim now is to compare the time evolution generated by the original quantum 
Hamiltonian H with the one generated by the approximate Hamiltonian Hqlt). For this 
we will follow the method devised in |CR97j for the case without spin. The presence of spin 
requires some modifications that, however, are modest when the spin quantum number s is 
fixed. But for the clarity of the presentation, and to prepare for the more involved situation 
to be dealt with in the second semi classical scenario, we will now present the argument in 
some detail. 

As stated in section |21 the Hamiltonian H generates a unitary and strongly continuous 
one-parameter group U{t,to), if its symbol satisfies the ellipticity condition ()2.5|) . When 
considering the limit ^ — *• and keeping s fixed this requirement need only be imposed on 
the principal symbol, i.e. we demand 

\Hoix,0 + i| > c (1 + \x\'Y'^^\l + . (4.11) 

Let now tjQ{t, to) be the corresponding unitary group generated by HQ{t). Using Duhamel's 
principle we may then express the difference between these unitary operators as 

U{t,to)-UQ{t,to) = ^ [ U{t,t'){H - HQ{t'))UQ{t',to) dt' . (4.12) 

J to 

Since we are interested in the difference (jH.lll . we have to consider the action of ()4.12|1 on 
the initial state </2^p) ® (f>n with to = 0. This requires an estimate of 

\\{H ~ HQ{t'))Mt')\\ , (4.13) 



14 



where ipQ^t') is the time dependent coherent state ()4.3|) . One can achieve this with the 
help of the following lemma, which is an immediate extension of a result given in jnR92|. 

Lemma 4.2. Let f,g e C°°(T*E'^) be symbols that satisfy the estimate \2.'J\) with M = 
and let F : T*E'^ ^ T*E'^ be a linear map with Hilbert- Schmidt norm ||F||hs- a,j3 & 
fi^'^ with k := \a\ = \P\ + 2 > 2 and introduce the symbol 

A{w) := (FwyfiFw) + h{FwYg{Fw) . 

Then for any real number k > there exist C > and G IN such that 

||op^[A]^,|| <Cn^/M llFllls sup \dlf{w)\ + \\Ft^^' sup \dlg{w)\ 

\ |7|<fc+iV |7|<fc-2+Af 

holds for any function ipfi(x) = fV^/^ipix/ \/h}j with if) G ^(R°') and < /i+ v^||-^||hs < 

We intend to apply this lemma to the difference ()4.13|) . with / corresponding to the 
Taylor remainder of Hq of order three and g to the Taylor remainder of C ■ d7rs(cr/2) of 
order one. But first we replace ()4.13|1 by 

l|f/Q(t',0)*(^-^Q(t'))f/Q(t',0)V^Q(0)|| (4.14) 

and invoke an appropriate Egorov theorem. Since the Hamiltonian generating f/Q(t,0) 
has a symbol that is composed of a scalar and quadratic principal part as well as a matrix 
valued subprincipal part, one can combine the techniques used in |BG00j and |Sch01j . This 
shows that 

W{t) ■.= UQ{Q,t){H - HQ{t))UQ{t,^) (4.15) 
is a Weyl operator with symbol 

W{t,w) = d*{z{t)) {H - HQ{t)){z - S^iit){w - z{t))) d{z{t)) . (4.16) 

Here d(^z{t)) is the representation ns{g{t)) of the solution to equation ()3.10|1 in which C{t) 
stands for C(^z{t)). Thus 

rf(z(t))0„ = e^^^W 

describes the transport of a spin-coherent state along the trajectory z{t). Since the principal 
part of the symbol H — Hqit) is scalar it is not affected by the conjugation with d[^z(t)) . 
In the subprincipal term this conjugation rotates the spin operator S = hd'iTs{a/2) to 
R(^g(t))S. Therefore, the spin part of the Egorov relation ()4.16p does not contribute to an 
estimate of ()4.14|1 in an essential way. 

If one now localises the symbol ()4.1(j|l in w with some smooth function that is compactly 
supported around z(t), leading to an error of size 0{h°°) when one applies W to a coherent 
state located at z{t), one can proceed to use Lemma 14.21 as in |(]R97j . This shows that 
there exists a constant K > such that 

\\{H - Hgrngm < K0'e{tf5{tr , (4.17) 
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where 

^(t) := maxjl, sup ||S'o,^(t')l|Hs| and := sup (l + (4.18) 

t'e[o,t] ' J t'e[o,t] 

depend on the classical trajectory z(t) = {q{t),p(t)). The constant m = ma,x{Mx:, M^} is 
related to M = {M^, M^) appearing in ()2.3p . We then obtain: 

Theorem 4.3. Let the conditions imposed on the Hamiltonian in section\^and the ellip- 
ticity condition ( |y/.ii| ) hold. Then the coherent state ipQ{t) defined in ^■!^ semiclassically 
approximates ip{t) = U(t, 0){'ffqp) ® <Pn) in the following sense, 

\m)-Mm<KVhte{tf . (4.19) 

The right-hand side vanishes in the combined limits h ^ and t ^ oo as long as t <^ 
Tz{h). The time scale Tz{h) depends on the linear stability of the trajectory z{t). If the 
latter possesses a positive and finite maximal Lyapunov exponent \raa.x{z), one has Tz{K) = 
g^^-^-^l log/i|. In the case of a trajectory on a (non-degenerate) KAM-torus this time scale 

isT,{h) = Ch~^l^. 

Proof. Conservation of energy, HQ[z{t)) = E, together with the ellipticity condition ()4.11|) 
implies that 6{t) is bounded from above by some constant depending on E. Thus the 
estimate ()4.17|) immediately yields ()4.19|) when used in ()4.12p . 

If z{t) is a trajectory with a positive, but finite, maximal Lyapunov exponent the 
dominant behaviour as t — > oo comes from the term 6{t)^. This is due to the relation 

Amax(^) = limsupy log||S'o,^(t)||HS , 

see (IA.1|) . which readily implies Tz{h) = g^;-^-^| log ^|. In the appendix we also discuss 
sufficient conditions under which finite maximal Lyapunov exponents occur. 

If z{t) is a trajectory on a KAM-torus one can introduce local action-angle variables 
(/, 0) in a neighbourhood of that torus such that in these canonical coordinates the fiow 
reads I{t) = I and (f){t) = + uj{I)t, see |Laz98j . One therefore finds 

\\SoAmls = ^d + f{I)t' , 

such that 9{t) ~ Kt as t — > oo, which finally yields Tz{h) = C h~^^^. In the degenerate 
case, where /(/) = 0, this changes to Tz{h) = C IT^/'^. □ 

In a next step we want to improve the semiclassical error in ()4.19|) to an arbitrary 
(half-integer) power of h. This requires higher order approximations that may be achieved 
as in |(yR97j by iterating Duhamel's principle ()4.12|1 . resulting in the Dyson expansion 



f/(t,0)-f/Q(t,0) = ^(in)-M .../ UQ{t,0)Wit,)...W{t,)dt,...dt, 

Jo Jtj^i (4.z(Jj 



+ RN{t]h) 
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with remainder term 



Jo JtM-l 



In order to estimate the contribution of the remainder when ()4.2n|l is apphed to the initial 
coherent state '(/'(O) = vfqp) ® 4'n we use the argument leading to ()4.17j) repeatedly. This 
yields the bound 

\\RN{t;h)^{Q)\\ < Kr, h^'/H'' 9{tf'' S{tr^ . (4.21) 

We then replace the symbol of each difference H — Hqitk) appearing in the sum in ()4.20|) 
by its Taylor expansion, 

!/|=3 ■ \u\=l 

(4.22) 

The integers Uk are chosen sufficiently large such that, after quantisation, the contribution 
of the remainder to an application of ()4.2()|1 to ip{0) can be absorbed in the error estimate 
fl4.21|) . Similar to the case without spin treated in |CR97j the quantisation of the main 
terms in ()4.22|) produces matrix valued differential operators Pkj(t) = op^[pkj(t)] with time 
dependent coefficients acting on the coherent state '^fqp-j ® 4>n- The symbols Pkjit){x,^) 
are polynomials in (x,^) of degree < k. Lemma f4. 21 finally leads to the following result: 

Theorem 4.4. Suppose that the quantum Hamiltonian H with symbol 112. ^} satisfies the 
conditions specified in section]^ and the ellipticity condition (|^.ii| ). Then for t > and 
any N e ¥1 there exists a state ipN{t) G L^(E"') ® C^'*^^ localised at (^q{t) , p{t) , n{t)) , that 
approximates the full time evolution ip(t) = U{t,0)[(p^^p^ ® 4>n) of a coherent state up to 
an error of order h^/'^ . More precisely, 

N-l / ,x j 

i=i ^ ^ 

The right-hand side vanishes in the combined limits h ^ and t ^ oo as long ast ^ Tz{h), 
where Tz{K) denotes the same time scale as in Theorem \4.^ 

Furthermore, ipN{t) arises from ^pfgp) <S> (pn through the application of certain (time 
dependent) differential operators Pkj(t) = op^[pkj(t)] of order < k, followed by the time 
evolution generated by Hqit), according to 

^^(t) =V^Q(t)+ f>Q(^,0)]5,,(t)^(0) . 

{fc,j)GAjv 

Here we have defined ■= {{kj) G IN x IN; 1 < k-2j < N-1, k > 3j, 1 <j< N-1}. 
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We remark that the matrix valued differential operators Pkj{t) do not increase the 
frequency set of a semiclassical distribution such as the initial state ffqp-j ® 4>n- This 
follows for the translational part from the respective statement without spin |R,ob87j . 
whereas the spin part is only acted upon by a matrix producing linear combinations of 
ipn- Moreover, according to Proposition 14.11 UQ(t,0) propagates the frequency set along 
the trajectory n(t)) so that both ipg^t) and ipwit) are semiclassically localised 

at (g(t),p(t),n(t)). 



4.2 Semiclassics with hs fixed 

We now consider the second semiclassical scenario in which both semiclassical parameters, 
h and s, are used. For this purpose we still represent the Hamiltonian H as a. matrix valued 
semiclassical Weyl operator. That way h appears as before, whereas the second parameter 
s G lN/2 controls the dimension of the space C^'^"'"^ on which the symbol operates as 
a linear map. As we will see, the parameter s enters relevant estimates through the 
expression hdiis{o'/2). To leading order this will produce factors of hs. Our desire to 
perform systematic semiclassical expansions therefore forces us to keep the combination 

S := hs 

fixed in the semiclassical limit. This means that from now on we consider ^ ^ and 
s — > oo with hs = S. 

An inspection of Proposition l4. H and its proof reveals that replacing hs by the constant S 
will shift the spin-action term g{t), which before was of subleading semiclassical order, to an 
additional contribution to the action Rq. This suggest that now the translational classical 
dynamics will be influenced by the spin, requiring a modified quadratic Hamiltonian. Not 
only that, revisiting the proof of Theorem 14. HI shows that we also have to estimate the 
application of spin operators to spin-coherent states in terms of s. This requires knowledge 
of the following: 

Lemma 4.5. For any X = x ■ cr/2 G su(2), n G and iV G IN there exist differential 
operators dH^ of degree 2j on C°°(S^) Cg> C^'*"*"^ and constants Cn > such that 



1 1^1 



< ^ • (4-23) 



i=o 

l(0) 



The leading order in this asymptotic expansion is determined by the constant Dn = 1, 

d7fs{x■a■/2)(j)r^ = sx■n(j)r^{l + 0{s'^)) . (4.24) 

Proof. We start with expressing a linear map L on the representation space C^^"*"^ in terms 
of Berezin's quantisation, 

L = {2s + l) [ P[L]{n)U{n) dn , (4.25) 
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where P[L\ denotes the upper (or P-) symbol of L, see e.g. jSim80| fPerSGj . Furthermore, 
dn is the normalised area measure on and T\.{n) stands for the projector onto the one- 
dimensional subspace in C^'^"^^ spanned by the coherent state vector (^n- In the present 
context the relevant linear maps are representation operators of Lie-algebra elements X = 
X ■ (t/2 E su(2). Their upper symbols are simple, 

P[d7r,(x ■ cr/2)](n) = {s + 1) x ■ n , 

see jSimSOl IPer86j , so that an application of such an operator to a coherent state reads 

dTTs{x-o-/2)(f)n= {2s + l){s+l) m ■ X {(pm, 4>n) 4>m dm . (4.26) 

The coherent states not being defined globally on is irrelevant to this expression since 
these states have been defined on a set of full measure. 

An asymptotic expansion of the integral ()4.2(jj) . as s — *• oo, can be achieved with the 
method of steepest descent. This is a variant of the stationary phase method, with a 
complex phase function, and is described in detail in |Hor90j . The first step consists in 
identifying the relevant phase factor, which in the present case is given by 

(0r.,0rr^)=e^^^"("^) with Im y.„(m) = - log (^^^^-^^ , (4.27) 

where ipn is independent of s, see |Per86j . Outside of a neighbourhood of m = —n the 
function limpn is finite and non-negative; it has a unique minimum at m = n. The real 
part of the phase cpn can be identified as the area of the spherical triangle with edges defined 
by the north pole, n and m. Hence m = n is the unique, non-degenerate stationary point 
of the phase. Up to an error of size ©(e"**) one can hence cut out a neighbourhood of 
m = —n from the integral ()4.26|) and use the representation ()4.27|) for Imyj^. The method 
of steepest descent then implies the existence of differential operators Dn^ of order 2j on 
C°°(S^) (g)C^''+^ and constants Cat > such that for any iV e IN the expansion (j03|l holds. 
The constant Dn^ fixing the leading order can be identified by choosing x = n, since 

d'Ks{n ■ Cr/2)(j)n = S(f)n ■ 

Comparing with ()4.23j) therefore yields Dn^ = 1, which implies ()4.24|1 . □ 

When constructing a quadratic Hamiltonian we now have to take into account that an 
application of a spin operator to a spin-coherent state contributes to the leading semiclas- 
sical order, as ()4.24p means 

We are therefore led to define a quadratic Hamiltonian Hgit) = op^[HQ{t)] with matrix 
valued Weyl symbol as follows, 

2 2 

\v\=^ \v\=\ (4.28) 

^hC{z{t)) ■d7r,(o-/2) . 
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Like in ()4.H) we have introduced a yet to be determined trajectory z(t) = (^q(t),p(t)) in 
T*R'^ with initial condition z{0) = z = {q,p), as well as a curve n{t) on with n(0) = n. 
This Hamiltonian, being quadratic in {Q, P) and linear in S, propagates an initial coherent 
state exactly: 

Proposition 4.6. The solution of the quadratic Schrddinger equation 

^f^^it) = HQ{t)^Q{t) With MO)=vip)®<l>r^ (4.29) 

is, up to an additional phase, again a coherent state, 

V,«(0=eK^-i'''))<;*!,^,,,®^„,,,. (4.30) 

Here [q(t),p(t),n{t)^ = ^l^{p,q,n) is the solution of Hamilton's equations of motion 
on T*E'^ X generated by the classical spin-orbit Hamiltonian 

H^x, ^, n) := Ho{x, + Sn ■ C{x, • (4.31) 

The phase of il)Q{t) is determined by 

Rsoit) = !\p{t')q{t') - H,{q(t'),p(t'))) dt' + SQ{t) , (4.32) 

^0 

which can be viewed as a total spin-orbit principal function, and by the Maslov phase cr{t). 
The latter derives from the time evolution 

B{t) = {dgq{t)B + dpq{t)){d,p{t)B + dpp{t))-' (4.33) 

of the complex symmetric d x d matrix i? G S^. 

Proof. The proof of this proposition parallels that of Proposition 14.11 however, a few 
modifications are necessary. One can again consider ()4.30|) as an ansatz and determine 
its ingredients by inserting it into ()4.29p . leading to equations analogous to ()4.8p . As 
opposed to ()4.9|1 the fact that now S* = is fixed shifts the term with g{t) from the 
last equation to the first one. Moreover, due to the modified definition of the quadratic 
Hamiltonian the principal symbol Hq is replaced by Hso in all places but one, yielding 

9 = qp- Ho + Sg 
-p + Bq = Hi^.^ + BHl^^^ 

The first two equations fix the translational part of the classical dynamics to be solutions 
of ()2.8|) with some n{t) and yield the spin-orbit principal function ()4.32|) . In the last two 
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equations Hso, which is evaluated at (^q(t),p(t),n(t)), can be viewed as a time dependent 
Hamiltonian, Hso{w,t) = Hso{w , n{t)) , for the translational degrees of freedom, with the 
time dependence introduced through n(t). These equations can be solved in the same 
manner as in the time independent case, yielding 

B{t) = S,oAt)[B] 
as in ()4.10|) . Here Sso^zit) is a solution of 

^5so,.(t) = JKim^t) SsoAt) (4.34) 

with S'so,z(0) = l2d', it hence yields fj4.33p . 

The classical spin motion n{t) so far has remained undetermined. Since the equation 
for the spin-coherent state is again ()4.7j) . it follows that n{t) must be a solution to the spin 
part of dZBl). □ 

In contrast to the previous case the classical dynamics that governs the time evolution 
of the coherent state ipQ{t) now is Hamiltonian on the product phase space T*E'^ X see 
(I2.8|l . This means that the spin precession is not merely following the translational motion, 
but there occurs a mutual influence of both dynamics. This effect is caused by the energy 
scales of the translational and the spin dynamics being comparable in the semiclassical 
limit, whereas when s is fixed the energy scale of the translational motion dominates. 

We now compare the time evolution generated by the full Hamiltonian H with the 
approximate dynamics following from the quadratic Hamiltonian Hqit) whose symbol is 
given in ()4.28|1 . As opposed to the situation discussed previously, see fl4.11|l . when keeping 
hs fixed the ellipticity condition has to be imposed on the full symbol of H, see fl2.5j) . which 
implies 



c (1 + \x\y^'^^\i + > II {H{x,o + ri > 

Here in the middle || ■ || denotes the operator norm on C^*^^, and on the right-hand side ip 
is any non-zero vector in C^'^"'"^. Choosing ip = + ^^(pn and using ()4.24|1 we then 

conclude that the spin-orbit Hamiltonian 1)4.311) is elliptic, in the sense that 

\HUx,^,n) + i| > C (1 + \x\'f^^\l + 

holds for all (x,^,n) e T*R'^ x S^. Therefore, we can again base our further investigation 
of the difference between the two quantum dynamics on the Duhamel relation ()4.12|) . This 
requires to estimate the analogue of ()4.13p . where in the present situation H — Hqit) is 
the Weyl quantisation of the symbol 



{H - HQ{t)){w) = J2 ^{Mn^i^m - Suit)) ■ C^^\zit)) {w - z(t)) 

\u\=l " 

+ H^o^ (t, w) + hC^^^ {t, w) ■ d7r,(o-/2) 



(4.35) 
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in which Hq and ' denote Taylor remainders of order three. Introducing an operator 
W{t) as in ()4.15|1 . the same type of an Egorov theorem as above apphes, leading to the 
symbol 

W{t,w) = d*{zit)) {H - HQ{t)) {z - S;lit){w - z{t))) d{zit)) (4.36) 

of W{t). We remark that z{t) being the projection of ^l^{z,n) to T*R'^ here requires the 
differential Sso,z{t) of with respect to z. The conjugation with d(^z{t)) has no effect 
on the scalar terms in ()4.35j) . whereas it rotates the spin operator to R(^g{t))hd7is{(T/2). 
Hence, for the application of ()4.3(j|l to a spin-coherent state (pn we can employ Lemma [4.51 
By also converting estimates with respect to s into ones with respect to h this yields to 
leading order 

{R{g{t))hdn,{a/2) - Sn{t)) <Pr. = S {R{g{t)) n - n(t))0„ + 0{s-^) = 0{h) . (4.37) 

Moreover, the complete asymptotic series in powers of provided by Lemma f4.5l results 
in a full asymptotic expansion of ()4.37|1 in powers of h. This observation now enables us 
to apply Lemma 14.21 in a completely analogous way to that used previously, yielding 

||(i7-i7Q(t))^Q(t)||<ir;z3/2^(t)3 5(t)-. 

Here the quantities 6{t) and S{t) are defined as in ()4.18|1 . however, now with the differential 
Sso,z(t), and z{t) as given in Proposition I4.(i[ 

The stability of the trajectory z(t) is encoded in the quantity 

Amax(^) = limsupylog||S'so,^(t)||HS • (4.38) 

Since z{t) is not the integral curve of a flow, rather than calling Xmaxi^) a Lyapunov 
exponent we refer to it as a stability exponent. This can, however, be bounded by the 
maximal Lyapunov exponent of the flow-line {z{t),n{t)) in T*R'^ x S^, see the appendix. 
Thus, in close analogy to Theorem 14.31 we finally obtain: 

Theorem 4.7. Let the conditions imposed on the Hamiltonian in section\^hold. Then the 
coherent state 'ipQ{t) defined in \4-3G^ semiclassically approximates %l){t) = U{t,0)[ip^^p-j ® 
(pn) in the following sense, 

\\m-Mt)\\<KVhte{t)' , 

when hs is kept fixed. The right-hand side vanishes in the combined limits h 0, s —>■ oo 
and t —>■ oo as long ast<^ Tz{h). The time scale Tz{h) depends on the linear stability of the 
trajectory z(t). If the latter possesses a positive and finite stability exponent Amax(-2), one 
has Tzih) = ^ I log h\ . In case z{t) is a projection to T*R'^ of a trajectory {z(t), Ti(t)) 

on a (non- degenerate) KAM-torus in T*R°' x this time scale is Tz{h) = C hr^^^ . 
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As in the previous case an improvement of the semiclassical error can be achieved with 
the Dyson expansion ()4.20|) . The present case, however, requires an additional estimate of 
the spin contribution in terms of s. Concerning the error term RN{t', h)ijj{0), the transla- 
tional part is dealt with by a repeated application of the argument leading to Theorem 14. 71 
For the spin part an inspection of the relations ()4.H5|1 and ()4.H(ij) reveals the necessity to 
estimate the successive apphcation of the operators 

A(t,) := C^''\z{tk)) ■ (R{gih))hd7rs{(T/2) - Sn{U 



to the spin-coherent state (pn- Representing these operators in the form ()4.25|1 . the result 
of their /-fold {I < j) application reads 

A{ti)...A{h)(j)r, = {2s + iy j ... jp[A{ti)]{mi)... P[A(tO] (m,) x 

g2 g2 (4.39) 

X n(m/) . . . Il{mi)(f)n dm; . . . dmi , 

with the lower symbols 

P[A(tfc)](mfc) = C^^^zitk)) ■ (s{R{g{tk))mk - n{tk)) + hR{g{h))m^ . (4.40) 



Starting with m/, the integral ()4.39|1 can be successively evaluated with the method of 
steepest descent similar to the proof of Lemma 14.51 The relation 



Ii{mi) . . . n(mi)0„ = 0rri,_i) 



then shows that the critical points of the phase are given by rrii = mi_i = ■ ■ ■ = mi = 
n. At these points, however, the lower symbols P[A(tfc)](mfc) are of order h, compare 
()4.40|) . The application of the method of steepest descent therefore yields in leading order 
a contribution 0{f^) = 0(s~'). Derivatives of total order n contribute terms of the order 
0{s~^fi~^) = 0(s~'), if n < Z, and of the order 0{s~^) otherwise. Altogether there hence 
exist differential operators V^'^'' of order < 2k on C°°((S^)') ® C^'^"''"'^ such that 

Mti) . - (l + E ^V^-\P[Ammi) . . . ^[A(ti)](mO0.)^^^..._^_ 

K=l 

(4.41) 

is of the order s for any K > I. The left-hand side of ()4.39|1 hence is of the order 

0(s~') = 0{h^), meaning that every factor A{tk) contributes a factor of h. We therefore 
finally obtain an estimate of the remainder term to the Dyson series given by 

\\RN{t- h)ij{{))\\ < KNff'H''e{tf''6{t)^'' . 

The main terms in the Dyson expansion are treated by replacing each factor of ()4.35p . 
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occurring at t = t^, with the Taylor expansions 

2 1 

^-(M7r,(o-/2) - Sn{tk)) ■ C^^^zih)) {w - z{h)y 
\u\=i^- 

+ £ ^ i^tk)) + hC'^"') (^(4)) d7r,(o-/2)) {w - z{tk)Y + nitk, w) , 

iy=3 

where again the integers are chosen sufficiently large. The contribution of the transla- 
tional degrees of freedom can be dealt with as in the previous semiclassical scenario, and 
the spin contribution follows from the expansion ()4.41|1 . Finally grouping together terms 
of corresponding orders in h, we arrive at a statement analogous to Theorem 14 .41 

Theorem 4.8. Suppose that the quantum Hamiltonian H with symbol \2.^) satisfies the 
conditions specified in section [B Then for t > and any G N there exists a state 
ipNit) G L2(E'^) C^*"*"^, localised at (g(t), n(t)) , that approximates the full time 
evolution tjj{t) = U {t, 0) (v^^p) ® <Pn) of a coherent state up to an error of order h^/"^ when 
hs is fixed. More precisely, 

\\m-Mt)\\ <cnJ2U] i^mr^^'' ■ 

j=i ^ / 

The right-hand side vanishes in the combined limits h ^ 0, s ^ oo and t ^ oo as long as 
t ^ Tz{K), where Tz{h) denotes the same time scale as in Theorem ^. 7| 

Furthermore, ipN{t) arises from ff^p-j ® (pn through the application of certain (time 
dependent) differential operators quKjif) = op^[Pfcj(^)] ®^k., 

(fc,Kj)eAiv 

where pkj (t) is a polynomial in {x, ^ ) of degree < k and is a differential operator of order 
< 2k on C°°(S^) C^'^"'"^. Here we have also defined 

An := {{k, K,j) e 1 < k + 2k - 2j < N - 1, k + 2k > 3j, 1 < J < N - 1} . 

The semiclassical localisation of ipN{t) here is different from the situation covered by 
Theorem 14.41 in that the operators act on (pn. But these are differential operators and 
hence do not increase the frequency set. This means that ipNit) is semiclassically localised 
at $so(g, j9, n) and in in this respect is not different from the classically propagated coherent 
state ipQ(t). 

5 Discussion 

In the previous section we analysed the semiclassical behaviour of coherent states in two 
different limits. In various places we saw that the difference between the two cases is 
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expressed in the way the classical translational and spin motion are coupled. Otherwise 
the final results agree to a large extent. This includes the mechanisms of semi classical 
localisation in the product phase space T*R'^ x S^. 

The problem of how the localisation of an initial coherent state develops with time 
can be made more explicit by using semiclassical phase-space lifts of the coherent states. 
At t = the state ip{0) = V'^p) ® (f>n is concentrated in a neighbourhood of the point 
{q,p, n) G T*R"' x S^. This concentration can be measured in terms of expectation values 
{ip{0), AiIj{0)) of operators A = op^[A] that are quantisations of well localised symbols 
A G C^(T*E'^) M2s+i(C). For simplicity we also assume that A is independent of h. 
At later times iplt) can in both semiclassical scenarios be approximated by an appropriate 
coherent state ipqit), such that 

m),M{t)) = {Mt),MQ{t)) + o{i) , t<^TM- (5.1) 

The expectation value on the right-hand side has a phase-space representation 



{^Q{t),A^i^Q{t)) = yy^^^^iy[<(^;/](u;)(0„(,),A(ti;)0„(,))c2.+i d^ . (5.2) 

A comparison with clearly reveals that the state ipqit) is concentrated at the point 
{q(t) , p{t) , n{t)) in the semiclassical limit as long as the quadratic form GB{t)/h is strictly 
positive definite. Either of the time evolutions ()4.6p and ()4.33p of B now imply jSchOlj 

so that the spreading of ipQit) in T*E°', see ()3.6|) . is bounded according to 



< -. r. < 



trG'B(t) ||G'B(t)||Hs ||G^b||hs 

Here Sz{t) denotes the differential of the appropriate fiow with respect to {x,C,)- If z{t) 
now is a trajectory with maximal Lyapunov (or stability) exponent Xmaxiz) > 0, the 
requirement for the state ipQ{t) to remain localised therefore is t <^ 2A ^ {z) I ^1 ■ This 
time scale is three times larger than Tz{h), which is the estimated time in 1)5.11) for the 
coherent state ipqit) to still well approximate the full time evolution ^p{t). 

Let us remark that the limitations in ()5.1)1 . to approximate the expectation value in 
terms of a coherent state, derive from estimating the difference ip{t) — ipqit) in L^-norm. 
But the error term on the right-hand side of 1)5.1)) measures this difference in a considerably 
weaker form so that one might expect it to vanish as — > and t ^ oo also for times 
Tz{h) < t 3Tz{h). In the case without spin Bouzouina and Robert jBR02j proved that 



this indeed holds, suggesting that the same is true in the present setting. 

Expectation values in coherent states such as ()5.1)) can also be used to obtain the leading 
semiclassical description of the propagation of observables. To see this let A, as above, be a 
bounded Weyl operator and denote its quantum time evolution by A{t) = U(t, 0)*A U{t, 0). 
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Here, however, we do not necessarily require the symbol to be compactly supported. The 
relations ()5.1|) and ()5.2|) then remain valid so that 



imJmm = / l_^K1?k«^) (0.w,a(^)0^,))c2.+i dw+o{i) . 



Since A{t) is bounded it may also be expressed as a Weyl operator, with symbol A(t) such 
that for t -C Tz{h) equation ()5.1|) can be rewritten as 



{27rhY 



{27rhy 



B{0) 

Bit) 

^[Vz(t) ]i^) {(pn{t),A{w)(f)n(t))€-2s+i dw = o(l) . 



The semiclassical localisation properties of the coherent states discussed above therefore 
imply that in leading order the symbol of the time evolved observable A{t) can be expressed 
in terms of the symbol of A transported along the classical flow {q{t),p{t), ^(t)), 

{(j)n,A{t){q,p)(f)n)c2s+i - A(g(t),p(t))0„(i))c2.+i = o(l) . 

The C^*"'"^-expectation values in spin-coherent states are lower (or Q-) symbols (see e.g. 



Sim80| IPer86j ) of the matrix valued functions A{t) and A, respectively. In terms of this 



mixed phase space representation of operators, employing Weyl calculus for the transla- 
tional part and Q-symbols for the spin part, this means that the quantum time evolution 
of observables follows the classical dynamics in leading semiclassical order. This statement 
represents a limited version of an Egorov theorem and again is valid for both semiclassical 
scenarios discussed in the preceding section, up to the time scale t ^ Tz{K). 



Acknowledgments 

A major part of this work has been performed when both authors stayed at the Mathe- 
matical Sciences Research Institute, Berkeley. We would like to thank the MSRI for its hos- 
pitality and for the support extended to us. Financial support by the Deutsche Forschungs- 
gemeinschaft under contract no. Ste 241/15-2 is gratefully acknowledged. R.G. was also 
supported through the Doktorandenstipendium D/02/47460 by Deutscher Akademischer 
Austauschdienst . 



Appendix: Linear stability of Hamiltonian flows 

The flows $0 ^^'i introduced in section |2l are both Hamiltonian flows on symplectic 
phase spaces. They are generated by smooth Hamiltonian functions H on 2n-dimensional 
smooth manifolds M with symplectic forms uj. In the first case the Hamiltonian is Ho{x,C,), 
defined on the phase space M = T*R'^ so that n = d and u = dx A d^. In the situation 



26 



of classical spin-orbit coupling the Hamiltonian Hso{x,C,,'n) is given on M = T*R'^ x S^. 
Thus, n = d + 1 and u = dx AdC,+ dn, where dn denotes the normalised area two-form on 
the sphere S^. In this appendix we want to recall the notion of Lyapunov exponents and 
give sufficient criteria of their existence in terms of properties of the Hamiltonian function. 

The linear stability of a flow $* is determined by properties of the differential D<I>*(q;) 
which is a linear map from the tangent space Tq,M to T<j,t(Q,)M. It, moreover, is a multi- 
plicative cocycle over the flow $*, i.e. D$*+*'(a;) = D$*'($*(a)) D$*(q;). If one introduces 
a euclidean scalar product in the tangent spaces, this gives rise to the adjoint D$*(a)*. 
Then D$*(a)* D$*(a) is a non-negative symmetric linear map on Tq,M whose eigenvalues 
we denote by 

^^^'\a)>■■■>^^?-\a)>0. 
The 2n Lyapunov exponents of the flow $* at a G M are now given by the expressions 

Xk{z) := limsup^i^log/if'^z) , 

t— >oo -^t 

if these are flnite. The largest Lyapunov exponent Ainax(a) provides a quantitative measure 
for the linear stability of $* since it measures the leading rate of local phase space expansion; 
it can be obtained from the relation 

Amax(a) = limsup;^logtr(D$*(a)*D$*(a)) . (A.l) 

t— >oo ■^t 

Hamiltonian flows leave the energy shells 

Qe ■■= {a G M; H{a) = E} 

invariant. If is a regular value of the Hamiltonian function H, the energy shell fi^; is a 
smooth submanifold of M of dimension 2n — 1. In such a case two Lyapunov exponents 
are always zero. They correspond to the direction of the flow and the direction transversal 
to the energy shell. Of the remaining 2n — 2 Lyapunov exponents half are non-negative (if 
they exist) and the rest of the Lyapunov spectrum is given by minus the flrst half. 

In general it is not known whether the Lyapunov exponents are flnite. If, however, an 
energy shell il^; is compact, one can introduce the normalised Liouville measure as a flow 
invariant probability measure on Qe- In this case one can apply Oseledec' multiplicative 
ergodic theorem to the restriction of $* to this energy shell |Use68j ; it guarantees that the 
Lyapunov exponents are flnite for almost all points ohQe with respect to Liouville measure. 
Moreover, if the flow is ergodic with respect to Liouville measure Xk{c() is constant on a 
set of full measure. Since we want to consider also non-compact energy shells we now give 
alternative sufficient criteria for the flniteness of Lyapunov spectra. 

Proposition A.l. Let H G C°°(M) he a Hamiltonian function such that the Hilbert- 
Schmidt norm ofD^H is bounded on the energy shell ^Iem that contains the point a G M. 
Then the Lyapunov exponents Ai(a), . . . , A2n(«) are finite. 
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Proof. Fix a G M and introduce canonical coordinates {q,p) & U G R" x R" in a neigh- 
bourhood of a. Then in this neighbourhood D^H is represented by the matrix H"{q,p) 
of second derivatives with respect to {q,p)- In these coordinates we denote the flow by 
^*{q,p)] its differential satisfies the equation 

jB^\q,p) = JH"{^\q,p)) D<l*(g,p) , D<l*(g,p)|i=o = Isn , (A.2) 

where J = ( _ ^q" ) • -^y integrating ()A.2jl and taking the Hilbert-Schmidt norm one 
obtains ^ 

l|D$*(g,p)||HS<2n+ / || Ji7"(|.^(g,p)) ||hs ||D$^(g,p)||HS ds . 

Jo 

For simplicity we here assume that for the points $^(a;), s G [0,t], one can use the same 
system of canonical coordinates. Gronwall's inequality then yields the estimate {t > 0) 

||Dl'*(g,p)||Hs < 2nexp|t sup \\J H" {^'{q,p))\\Hs] < 2ne^* , 

with some constant C > 0. The last line follows from the boundedness of D^H on QE,a- 
Since on the other hand 

m\q,p)\\HS = \j ^^'l\a) + ■ ■ ■ + ^^f^\a) , 

the bound 

liog/ir»<^ 

for the maximal eigenvalue ^™^^{a) follows. This finally implies the assertion. □ 

An application of this Proposition to the two flows $q (defined on M = T*R°') and 
(defined on M = T*E'^ x S^) immediately yields 

Corollary A.2. // the norm of Hq is bounded on Qe,{x,^) C T*R'^, the 2d Lyapunov 
exponents Xo,k{x,^) of the flow $q are finite. If, in addition, the derivatives C'^''\x\^') of 
order |z/| < 2 are hounded for all {x',^',n') G flE,{x,(,n) C. T*E'^ x S^, the 2d + 2 Lyapunov 
exponents \so,k{x, ^, n) of the flow are also finite. 

In the second semiclassical scenario, however, rather than the Lyapunov exponent 
Ko,k{<l,P,'>^) of a point {q,p,n) G T*E'^ x the stability exponent ()4.38|) of the pro- 
jection to T*R'' entered Theorem 14.71 Revisiting the proof of Proposition I A. II shows that 
in view of ()4.H4j) such a stability exponent is finite under the same conditions as stated in 
Corollarv IA.2I for Aso,a:- Moreover, a simple estimate yields the bound 

Amax ^ Ago, max • 
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